Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN W^X style file vl.4) 



The AGN contribution to deep submillimetre surveys and the 
far-infrared background 

O. Almaini, 1 A. Lawrence, 1 B.J. Boyle 2 

1 Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ 
2 Anglo-Australian Observatory, PO Box 296, Epping, NSW2121, Australia 



On 

U 

< 
in 



> 
oo 

m 
O 

On 
On 

Or 

6 



Accepted 1999 March 24. Received 1999 March 15; in original form 1998 September 2 



ABSTRACT 

A great deal of interest has been generated recently by the results of deep submillimetre 
surveys, which in principle allow an unobscured view of dust-enshrouded star formation at 
high redshift. The extragalactic far-infrared and submillimetre backgrounds have also been 
detected, providing further constraints on the history of star formation. In this paper we esti- 
mate the fraction of these backgrounds and source counts that could be explained by AGN. 
The relative fractions of obscured and unobscured objects are constrained by the requirement 
that they fit the spectrum of the cosmic X-ray background. On the assumption that the spec- 
tral energy distributions of high redshift AGN are similar to those observed locally, we find 
that one can explain 10 — 20 per cent of the 850/im SCUBA sources at lmJy and a similar 
fraction of the far-infrared/submillimetre background. The exact contribution depends on the 
assumed cosmology and the space density of AGN at high redshift (z > 3), but we conclude 
that active nuclei will be present in a significant (though not dominant) fraction of the faint 
SCUBA sources. This fraction could be significantly higher if a large population of AGN are 
highly obscured (Compton-thick) at X-ray wavelengths. 
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1 INTRODUCTION 



Deep submillimetre observations offer the potential to revolutionise 
our understanding of the high redshift Universe. Longward of 
100/^m, both starburst galaxies and AGN show a very steep de- 
cline in their continuum emission, which leads to a large negative 
K-correction as objects are observed with increasing redshift. This 
effectively overcomes the 'inverse square law' to pick out the most 
luminous objects in the Universe to very high redshift (Blain & 
Longair 1993). Since the commissioning of the SCUBA array at 
the James Clerk Maxwell Telescope a number of groups have an- 
nounced the results from deep submillimetre surveys, all of which 
find a high surface density of sources at 850/im (Smail et al. 1997, 
Hughes et al 1998, Barger et al. 1998, Eales et al. 1998, Blain et 
al. 1999b). The implication is the existence of a large population 
of hitherto undetected dust enshrouded galaxies. In particular, the 
implied star-formation rate at high redshift (z > 2) could be sig- 
nificantly higher than that deduced from uncorrected optical-UV 
observations (Hughes et al. 1998). The recent detections of the far- 
infrared/submillimetre background by the DIRBE and FIRAS ex- 
periments (Puget et al. 1996; Fixsen et al. 1998; Hauser et al. 1998) 
provide further constraints, representing the integrated far-infrared 
emission over the entire history of the Universe (Dwek et al. 1998). 
Since most of this background has now been resolved into discrete 



sources by SCUBA, the implication is that most high redshift star- 
forming activity occurred in rare, exceptionally luminous systems. 

In this paper we investigate the possible contribution from 
AGN to both the far-infrared/submillimetre background and the 
SCUBA source counts at 850/^m. The deep SCUBA sources in par- 
ticular are believed to be the high redshift equivalents to local Ultra- 
Luminous Infrared Galaxies (ULIRGs; Sanders & Mirabel 1996) 
and hence a large AGN fraction may not be a surprise. Our strategy 
is to predict faint submillimetre counts based on our knowledge of 
the AGN luminosity function and its evolution, together with an as- 
sumed spectral energy distribution (SED). We take account of the 
likely population of obscured AGN, using models which reproduce 
the spectrum of the hard X-ray background (XRB). 



2 OBSCURED AGN AND THE X-RAY BACKGROUND 

Locally it is clear that a large fraction of AGN are obscured by gas 
and dust. Obscured (e.g. narrow-line) AGN are believed to outnum- 
ber unobscured (e.g. broad-line) objects by a factor of ~ a few, with 
significant uncertainties depending on the selection techniques and 
the assumptions made (Lawrence 1991; Huchra & Burg 1992; Os- 
terbrock & Martel 1993). It is now important to understand whether 
this ratio continues to high redshifts and luminosities. Radio emis- 
sion is unaffected by obscuration, and among radio-loud objects 
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at least there is clear evidence that narrow-line radio galaxies are 
common to very high powers and redshifts. Furthermore, recent 
deep X-ray surveys using ROSAT, ASCA and Beppo-SAX have de- 
tected substantial numbers of 'narrow-line X-ray galaxies', many of 
which show very clear evidence for obscured AGN activity (Boyle 
et al. 1995; Almaini et al. 1995; Ohta et al. 1996; Iwasawa et al. 
1997; Boyle et al. 1998; Schmidt et al. 1998, Fiore et al. 1999). It 
seems likely that this is the 'tip of the iceberg' of a large population 
of obscured AGN, full confirmation of which will soon be possible 
with the next generation of X-ray satellites (e.g. AXAF and XMM). 

There is increasing evidence that these obscured AGN are re- 
sponsible for the production of the hard XRB. At soft X-ray ener- 
gies (below 2keV) almost all of the XRB has been resolved into 
discrete sources, most of which turn out to be broad-line QSOs 
(Shanks et al. 1991), but the X-ray spectra of ordinary QSOs are 
too steep to explain the XRB at higher energies. The energy density 
of the XRB actually peaks at ~ 30keV, where ordinary broad-line 
QSOs can account for only ~ 20 per cent (Fabian et al. 1998). Ob- 
scured AGN provide a very natural explanation, since photoelectric 
absorption allows only the hard X-rays to penetrate (Setti & Woltjer 
1989). Models have been developed which provide very good fits to 
the XRB spectrum, X-ray source counts and the local distribution 
of absorbing columns (Comastri et al. 1995; Miyaji et al. 1998). 
The implication is that most of the energy density generated by ac- 
cretion in the Universe takes place in obscured AGN. As outlined 
by Fabian & Iwasawa (1999) this hidden population could explain 
the apparent discrepancy between predicted present day black hole 
densities and the observations of Magorrian et al. (1998). 



3 THE AGN LUMINOSITY FUNCTION 
3.1 Type 1 AGN 

Our predictions are based on the X-ray Luminosity Function (XLF) 
of Boyle et al (1994), which we take to represent the baseline pop- 
ulation of unobscured (Type 1) AGN. Since the optical/UV spectra 
of these AGN all show broad emission lines, and the X-ray spectra 
show no evidence for photoelectric absorption (Almaini et al. 1996) 
we can safely assume that this sample includes only relatively un- 
obscured objects (e.g. Nh < 10 22 cm _1 ). This luminosity function 
was modelled with a broken power law of the form: 



9 X Ll 44 



v x T — 



Lx < L* x 
Lx > L* x 



(1) 



where 71 and 72 represent the faint and bright end slopes of the 
XLF respectively and L44 is the 0.3—3.5 keV luminosity expressed 
in units of 10 44 erg s _1 . A good fit to the evolution of this luminos- 
ity function was obtained with a Pure Luminosity Evolution (PLE) 
model out to a maximum redshift z max : 



L*x{z) 



L x (0)(l + z) k z<z„ 



Lx {Zmax) 



(2) 



z > z„ 



The best fitting parameters were z ma x = 1.6(1.79), k — 
3.25(3.34), 71 = 1.36(1.53), 72 = 3.37(3.38), L* x = 
10 43.57(43.70\ and &* x = 1.59(0.63) x 10~ 6 Mpc~ 3 for a cos- 
mology with q () — 0.5(0.0) and Hq = 50 km s _1 Mpc _1 . Further 
details can be found in Boyle et al. (1994) (models S & T). 



3.2 Correction for obscured AGN 

To constrain the relative fractions of obscured and unobscured ob- 
jects we use the Comastri et al. (1995) population synthesis model 
for the XRB. In this model a population of obscured AGN with 
a range of column densities is added to the unobscured popula- 
tion in order to obtain good fits to the XRB spectrum and the X- 
ray source count distribution. Excellent fits to the hard XRB were 
obtained by assuming a simple model in which the space density 
of obscured objects and the distribution in obscuring column den- 
sities are free parameters. The best fit was obtained with a ratio 
of obscured to unobscured AGN (Nh < 10 22 cm _1 ) in the range 
2.4 — 3.7, in very good agreement with local observations (see Sec- 
tion 2). Although the affects of the photoelectric absorption on the 
X-ray source counts are complicated, in the submillimetre regime 
the gas and dust is transparent. We therefore simply scale the un- 
obscured QSO luminosity function as follows: 



^obscured — 3 X ^unobscured 



3.3 The high redshift evolution 



(3) 



At the highest redshifts (z > z max ) the X-ray luminosity function 
of Boyle et al. (1994) is formally consistent with being constant, 
but beyond 2 = 3 there are very few QSOs in this survey. Optical 
surveys for bright QSOs find evidence for an exponential decline 
in the QSO space density towards high z, parameterised as follows 
by Schmidt, Schneider & Gunn (1995): 

$(z) = $(2.7)e (2 ' 7 - z) 



z > 2.7 



(4) 



Similar results were found by Hewett et al. (1993). The re- 
ality of this exponential decline is unclear, however. We note that 
these optically derived surveys may be seriously underestimating 
the high redshift QSO population because of intervening absorp- 
tion along the line of sight. This is discussed further in Section 5. 
We therefore consider two models for the high redshift evolution. 
For Model A we use the exponential decline given by Equation 4, 
while for Model B we consider the possibility that the space density 
of QSOs is constant to z = 5 (with an exponential decline there- 
after), in good agreement with the findings of the Ultra-Deep X-ray 
Survey by Hasinger (1998). 



4 THE THERMAL FAR-INFRARED SPECTRAL 
TEMPLATE 

In radio-quiet AGN there is significant evidence that the far- 
infrared/submillimetre emission is dominated by thermal re- 
radiation from dust (Carleton et al. 1987, Hughes et al. 1993). To 
model this spectrum, we use the best fitting spectral energy dis- 
tribution obtained by Hughes, Davies & Ward (1999), who have 
obtained the largest collection of far-infrared and submillimetre ob- 
servations of local AGN. They find that the 50 — 1300pim emission 
is consistent with thermal re-radiation from dust at a temperature of 
35 — 40K. This can be modelled by an isothermal grey-body curve: 



fv oc 



,3+/3 



A > 50/im 



(5) 



exp(/iz//fcT) - 1 

in which T is the temperature of the dust and (i is the emissiv- 
ity index, derived by assuming that the grey-body is transparent to 
its own emission. Best fitting values were found to be /3 = 1.7±0.3 
and T = 37 ± 5K (see Hughes, Davies & Ward 1999 for full de- 
tails). We note that strikingly similar temperatures have recently 
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been reported for high redshift AGN by Benford et al. (1998). At 
mid-infrared wavelengths they approximate the spectrum with a 
power law of the form fv tx v~ 1A . 

In order to normalise these SEDs to our X-ray luminosity 
function, we use the mean far-infrared to X-ray relationship de- 
termined for a large sample of broad-line, radio-quiet quasars by 
Green, Anderson & Ward (1992). Similar ratios were found by Car- 
leton etal. (1987). 

/lOOpm = 3.2 X 10 J /lfe e v (6) 

We note however that locally there is a large scatter in this 
X-ray to far-infrared ratio. To allow for this we adopt a Gaussian 
distribution in the log of the flux ratio given above with a dispersion 
a = 0.19, as observed by Green, Anderson & Ward (1992). 

A major source of uncertainty in our method is clearly the ex- 
trapolation of these local SEDs to AGN at high redshift. Accurate 
far-infrared and submillimetre photometry of high redshift, X-ray 
selected QSOs would help to resolve this issue, which will cer- 
tainly be possible with SCUBA and future submillimetre detectors. 
If these AGN contain more isotropically distributed dust they may 
well turn out to be brighter at these wavelengths, which could boost 
the source predictions significantly. 



5 THE SUBMILLIMETRE COUNT PREDICTIONS 

Using the SED defined above, combined with the evolving X-ray 
luminosity function from Section 3, one can predict the AGN sub- 
millimetre counts at 850/im. The results are shown in Figure 1 
for two different cosmologies, where we compare with the total 
SCUBA detections from recent deep surveys. We note that a more 
recent model for the XRB has recently been produced by Miyaji 
et al. (1998), based on the Ultra-Deep Lockman Hole survey of 
Hasinger (1998). In this survey they find that a Luminosity Depen- 
dent Density Evolution (LDDE) model provides a better fit to the 
evolution of the AGN luminosity function then a PLE model. Using 
model parameters provided by Miyaji (private communication) we 
find submillimetre counts predictions which are approximately 40 
per cent higher than shown in Figure 1 at the flux densities of the 
SCUBA surveys. 

In using these luminosity functions, we integrate over the lu- 
minosity range 10 42 < Lx < 10 48 erg s _1 . Integrating to higher 
luminosities make no significant difference, due to the very low 
space density of such sources. There is, of course, significant evi- 
dence that AGN activity continues to fainter luminosities (e.g. Ho, 
Filippenko & Sargent 1997) but this also makes no difference to 
our predictions. At an 850^im limit of lmjy there is no significant 
contribution from QSOs below a luminosity of L x ~ 10 43 erg s~ 1 . 

A major source of uncertainty in determining these predictions 
is the space density of AGN at high redshift (z > 3). For Model 
A we adopt the conservative assumption that $ declines exponen- 
tially above redshift z — 2.7, in the same manner as determined 
for bright optical QSOs (Schmidt, Schneider & Gunn 1995). To il- 
lustrate the importance of the high redshift space density, we also 
plot the extreme case where <E> does not decline beyond z — 2.7, 
but instead remains constant out to a maximum redshift of z — 5 
(Model B). The exact behaviour in the QSO space density at high 
redshift makes little difference to the XRB but due to the negative 
K-correction the predictions in the submillimetre waveband are al- 
tered considerably. As discussed in Fall & Pei (1993), because of 
dust along the line of sight in damped Lya absorption systems, ex- 
isting optical QSO surveys may be seriously incomplete at high 



redshift, with possibly up to 90 per cent missing at z = 4. Radio 
emission is unaffected by such absorption, and the high redshift de- 
cline seen in the radio-loud population would seem to argue against 
any serious incompleteness (Shaver et al. 1996, Dunlop 1997). We 
note, however, that recent results from the deepest X-ray surveys 
seem to suggest that the space density of QSOs remains constant 
out to z ~ 5 (Hasinger 1998). We conclude that the exact form of 
high redshift evolution remains uncertain, but the reality is likely 
to lie somewhere between Models A and B. It is worth noting that 
the n(z) distribution of submillimetre selected quasars will in prin- 
ciple allow a very sensitive discriminant between these competing 
models (see also Blain et al. 1999a). 

It is clear from Figure 1 that a sizeable fraction of the existing 
submillimetre sources can be explained by AGN. Taking the lmJy 
detection by Blain et al. (1999b), for example, the q = 0.0 Models 
A and B account for 9 and 19 per cent of the sources respectively. 
The qo = 0.5 models give very similar source counts at brighter 
fluxes but flatten considerably at the faint end, accounting for only 
7 and 13 per cent of the counts at lmjy. The uncertainties in our 
method could account for a factor of ~ a few in either direction, but 
nevertheless we conclude that a significant (though not dominant) 
fraction of the recent SCUBA detections will contain an AGN. As 
discussed in Section 7, however, these predictions are almost cer- 
tainly lower limits since they do not include Compton-thick AGN. 



6 AGN AND THE FAR-INFRARED/SUBMILLIMETRE 
BACKGROUND 

As discussed in Section 2, there is now growing evidence that AGN 
could account for the majority of the X-ray background. Recent 
analysis of COBE FIRAS data has also led to the first detections of 
an extragalactic component to the far-infrared/submillimetre back- 
ground (100 - 1000/xm; Puget et al. 1996, Fixsen et al. 1998). 
At present there are considerable uncertainties involved in the sub- 
traction of galactic emission, zodiacal light and the cosmological 
microwave background, but by a variety of independent techniques 
Fixsen et al. (1998) find a mean spectrum for the far-infrared back- 
ground which can be parameterised by the following function: 

I v = 1.3 ± 0.4 x 10 _5 (^ o ) (, ' 64±0 ' 12 P„(18.5 ± 1.2JT) ( 7 ) 

where v = 100cm -1 and P v is the Planck function. The 
nature of this radiation is still unclear, but at 850/im most of 
the background has now been resolved into discrete sources by 
SCUBA (e.g. see Blain et al. 1999b). Using the models out- 
lined above we can explicitly calculate the potential contribution 
from galaxies containing bright AGN. Extrapolating down to an 
850/im flux limit of O.OlmJy the resulting background contribu- 
tions are shown in Figure 2. With a conservative estimate for the 
space density of AGN, one can explain ~ 10 per cent of the 
far-infrared/submillimetre background. With a higher space den- 
sity of AGN (Model B) the contribution increases, particularly to- 
wards longer wavelengths, with ~ 20 per cent of the background 
at 850/im. 



7 COMPTON-THICK AGN 

In the arguments above, we constrain the relative contributions of 
obscured and unobscured AGN using population synthesis models 
which reproduce the spectrum of the XRB. Above a critical column 
density of ~ iq 24 -25 atom C mT 2 , however, the obscuring medium 
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Figure 1. Predicted AGN number counts compared with recent SCUBA measurements of the total submillimetre source density at 850/^m. The unfilled 
circle is taken from an estimate of the source density at lmjy based on a fluctuation analysis (Hughes et al. 1998). In Model A (dashed lines) we adopt the 
conservative assumption that the AGN space density, 4>, declines exponentially beyond z = 2.7. Model B (dotted lines) assumes a constant $ out to z = 5. 



can become optically thick to X-rays by Compton scattering, lead- 
ing to almost total absorption. Recent BeppoSAX hard X-ray ob- 
servations of a well defined sample of Seyfert 2 galaxies have re- 
vealed that the absorbing columns are in general very much higher 
than previously thought, and that most appear to be Compton-thick 
(Maiolino et al 1998). The contribution of such objects to the XRB 
will be negligible, and hence the space density of AGN required 
to fit the X-ray background will always be a lower limit. Hence the 
predictions shown in Figures 1 & 2 could be substantially higher, by 
~ 50 per cent if we adopt the distribution found by Maiolino et al. 
(1998). We note that Compton-thick AGN appear to be very com- 
mon among local AGN-dominated ULIRGs (Sanders & Mirabel 
1996) which show very high central gas densities. 



8 CONCLUSIONS 

There is growing evidence that obscured AGN are responsible for 
most of the cosmic XRB. Using models for the XRB to constrain 
the relative fractions of obscured and unobscured AGN, combined 
with recent measurements of their far-infrared and submillimetre 
SEDs, we estimate the contribution of AGN to the far-infrared and 
submillimetre background. We also explicitly calculate the pre- 
dicted source counts at 850/im, where a lot of excitement has been 
generated recently from the results of deep SCUBA surveys. 

Several lines of argument lead to the conclusion that a signif- 
icant fraction of the SCUBA detections will contain bright AGN. 
The very high far-infrared luminosities implied for these sources 
suggest that they are high redshift equivalents to local ULIRGs. 
At the luminosities of the recently discovered SCUBA sources 
(L ~ 10 12 L o ) we note that 20 - 30 per cent of local ULIRGs 
show clear evidence for an AGN (e.g. Seyfert 1 or Seyfert 2), with 
a further similar fraction showing more ambiguous LINER activity 
(Sanders & Mirabel 1996). Our explicit calculations suggest that 



the AGN contribution in recent SCUBA surveys is likely to be at 
least 10 — 20 per cent, with the exact contribution depending on the 
assumed cosmology and the space density of AGN at high redshift. 
This estimate is constrained by the requirement that we fit the spec- 
trum of X-ray background, and hence the fraction could be much 
higher if a large population of AGN are Compton-thick, as recently 
observed by BeppoSAX for local Seyferts (Maiolino et al. 1998). 
Such AGN will may not be revealed even by AXAF or XMM. 

If the dust in these AGN is heated by accretion processes, 
a significant AGN contribution may help to explain the apparent 
contradiction implied by the recent SCUBA surveys, namely that 
such massive star-formation at early epochs could very easily over- 
predict the abundance and metallicity of low-mass stars in the local 
Universe (Blain et al. 1999a). In many other ways a large AGN frac- 
tion is not a surprise. Of the very few high redshift submillimetre 
sources with good optical spectra, at least 3 show clear evidence for 
the presence of an AGN (e.g. Ivison et al. 1998, Irwin et al. 1998, 
Serjeant et al. 1998). 

The most puzzling question is the origin of the far-infrared 
emission in AGN themselves, even if they do explain a large 
fraction of the SCUBA sources. At these wavelengths the rel- 
ative roles of AGN and starburst activity remain controversial, 
both in IR selected ULIRGs (Sanders & Mirabel 1996) and in 
quasars (Lawrence 1997). Is the dust heated by the AGN or a nu- 
clear starburst? Both energy sources are clearly present in a large 
fraction of these objects. At mid-infrared wavelengths the AGN 
probably dominates but at 100/im the situation remains unclear 
(Sander & Mirabel 1996). Hence the emission seen in the far- 
infrared/submillimetre background, and from the recently detected 
SCUBA sources, could be entirely due to star forming activity, al- 
beit with the interesting implication that a large fraction of the star 
formation in the early Universe took place in the cores of galaxies 
containing active quasars. 
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Figure 2. Predicted AGN contribution to the far-infrared/submillimetre background. The dashed and dotted lines are based on the Models A and B respectively. 
The solid black line represents an analytical fit to the background spectrum obtained by Fixsen et al. (1998). 
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